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We consider the problem of improving noisy quantum measurements by suitable preprocessing 
strategies making many noisy detectors equivalent to a single ideal detector. For observables per- 
taining to finite-dimensional systems (e.g. qubits or spins) we consider preprocessing strategies 
that are reminiscent of quantum error correction procedures and allows one to perfectly measure an 
observable on a single quantum system for increasing number of inefficient detectors. For measure- 
ments of observables with unbounded spectrum (e.g. photon number, homodyne and heterodyne 
detection), the purification of noisy quantum measurements can be achieved by preamplification as 
suggested by H. P. Yuen [l[. 

PACS numbers: 



I. INTRODUCTION 

In many situations it is necessary to measure an observable in the presence of noise, e.g. when transmitting a 
quantum state through a noisy quantum channel that degrades it exponentially versus distance, corresponding to a 
degradation of the measurement. 

A number of figures of merit can be used to characterize the noise of non-ideal measurements. An example of 
such figures of merit is the variance of the outcomes distribution. An extensive analysis of the variance affecting 
quantum measurements has been done for example in Q. In a communication scenario, a relevant figure of merit is 
represented by the mutual information between the measurement outcomes and the input alphabet encoded on an 
ensemble of states. The problem of how much classical information can be extracted from a quantum system has 
been first deeply discussed by Holevo Q, who provided bounds on the accessible information, and then revisited in 
the framework of quantum information by Schumacher et al. [J] . A further figure of merit is the average probability 
of correctly distinguishing input states picked up from a given ensemble. This is one of the first problems faced by 
quantum estimation theory, and has been addressed extensively in the literature [H, [H-@] . Finally, another example of 
figure of merit is a suitable distance between the noisy and the ideal outcomes probability for fixed input states. 

In this paper, we consider the situation where N identical preparations of the state p g belonging to some ensemble 
S = {(p g , p g )} are given. We are allowed to use M non-ideal detectors, with M > N. Each detector is described by a 
Probability Operator- Valued Measure (POVM), namely a set of positive operators {P/}, which provides a resolution 
of the identity, i.e. JV P{ = P Each POVM element P- is the noisy version of an ideal POVM element Pi. A generic 
quantum channel 1Z is allowed to act on the N identical copies of the state p g before the M noisy POVM {P/} are 
measured, and a generic classical post-processing can be done on the outcomes of such measurements. Such a scheme 
of "purification" of noisy measurements is depicted in Fig. [TJ We address the problem of optimizing the quantum 




FIG. 1: Purification scheme for noisy quantum measurements. 

channel 1Z in order to reduce the effect of noise affecting the POVMs {P/}. We approach the problem through the 
minimization of the variance of the maximum likelihood estimator for the parameter g and through the maximization 
of the mutual information between g and the measurement outcomes. 
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Notice the analogy between quantum error correction schemes @, as depicted in Fig. [2J and the purification of 
measurements. For error correction, the message is first encoded by gate 1Z into one of the carefully chosen codewords, 
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FIG. 2: Scheme for quantum error correction. 



which is then (possibly) corrupted by the noisy communication channel £ . Finally, in gate T> some set of commuting 
Hcrmitian operators are measured over the corruption, the syndrome is used to perform error correction, and finally 
the recovered codeword is decoded into the original message. For purification of measurements, we are allowed to 
encode the N identical copies of input state p g through the channel 1Z, in a way similar to quantum error correction. 
The aim of such encoding is very different, since after that we are forced to perform M measurements with the same 
noisy POVM {P/}, which provide us just classical outcomes to be classically post-processed. The limitation of the 
measurement purification versus error correction is that the decoding T> is restricted to classical outcomes only. The 
problem we are considering is also similar to the problem solved by entanglement purification protocols Q, since we 
are generally trying to recast the use of a number of noisy measurements to an effective use of a smaller number with 
less noise. 

The paper is organized as follows. In Sect. |TT] we specify the general problem to a qubit with isotropic noise, 
and then we face the optimization considering different figures of merit: in Sect. Mil we show how to minimize the 
measurement noise, while in Sect. IIVI we maximize the mutual information between the parameter describing the 
state and the outcomes of the POVMs. In Sect. |V]and Sect. I VII we consider observables with unbounded spectrum, 
for which the concept of amplification applies, and we review the scheme of H. P. Yuen [l[ for purifying photodetectors 
(Sect. [V}, homodyne and heterodyne detectors (Sect. IVI[) . Finally, Sect. IVIII is devoted to conclusions. 



II. PURIFICATION OF QUBIT MEASUREMENTS 



Let us specify the general problem we are considering. We are provided with N identical copies of the input state 
p g of dimension d. In what follows we will always suppose that the elements of the POVMs {Pi} and |P/| are d, 
which has been proved to be the optimal choice for d = 2 [10| . when the mutual information is optimized 30]. We 
suppose that each noisy element P[ is obtained acting with the same channel £ on the corresponding element Pi of 
the ideal POVM 

where £ v denotes the Heisenberg-picture version of the channel £. Eq. (fTJ) shows that the ideal POVM {Pi} is 
"cleaner" that the noisy POVM {P/} in the sense of the partial ordering introduced in [l2j], as depicted in Fig. [3] 




FIG. 3: Noisy POVM element. 



We consider a qubit (so d = 2) parametrized as 

Pa,l = ( £ ! t a ) ■ (2) 

We are interested in the observable a z , and we suppose to have at our disposal M noisy POVM {P/} of a z> i.e. 
Pi = We assume a simple kind of noise acting on each POVM, i.e. the isotropic noise 

£ v {Pi) = aPi + pI, (3) 

so that P/ = a\i)(i\+/3I. 

We suppose to have N = 1 qubit state and consider as a purification channel 1Z the orthogonal cloning C, with 
respect to the basis of eigenstates of the observable a z 

C(p)= $>|p|*>|*><i|®". (4) 

i=0,l 
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The conditional probability p(i\a, b) of obtaining outcomes i — {ii, . . . , im} given the state parametrized by a, b 
does not depend on 6, and can be explicitly written as 

pflo) = Tr[C(p)S v (Pi)® M } . (5) 

We substitute Eq. @ and Eq. © into Eq. © to obtain 

p(i\a) = Tr[(a|0)(0|« M + (1 - a)\l){l\* M ) ®f =1 {a^)^\ + (31)} (6) 

We observe that the probability p(i\a) depends only on the number of outcomes O's and l's in the measurement (i.e., 
not on their position). Upon defining such integers as M$ and Mi = M — Mo, we obtain 

p(Mi\a) = [a(a + P) M °j3 M > + (1 - a)(a + fi) M ^ M °] . (7) 

For the normalization condition of the POVM in Eq. ([3]) one has a = 1 — 2/3, so < j3 < |, and hence 

p(M x |a) = W [a ((1 - - (1 - /3) M ^ M ») + (1 - Z?)^/?^] . (8) 

One can easily check the normalization of this probability, i.e. ^2m 1=0 p(Mi\o,) = 1. In the case of ideal measure- 
ments for which /3 = 0, the non-null probabilities are obtained just for Mi = and for Mi = M, namely 

p(Mi = 0\a) = a, p(Mi = M\a) = 1 - a, (9) 

whereas in the completely isotropic case (i.e. (3 = -|) the probability p(Mi\a) = {m^} (f) is independent of a, 
namely no information can be obtained about the state. 

Notice that also the coherent channel, widely used in encoding schemes for quantum error correction as (l3j . 

c'( P )= £ (i\ P \j)\i){j\® M . (io) 

ij=0,l 

leads to the same probability distribution Eq. ©, since P[ are diagonal on the a z basis. 

III. MINIMIZATION OF MEASUREMENT NOISE 

We show how to apply the ML criterion to obtain the optimal estimator for the expectation value of a z , by means 
of our measurement purification scheme. Our aim is to show an improvement of estimation in terms of variance by 
increasing the uses of the POVM. In the following log will denote the logarithm to the base 2. 

The ML criterion provides the following estimator for a in the state Eq. © 

Qm l = argmax — L(a\Mi) (11) 
a n 

where n is the number of (joint) outcomes (runs of the purification scheme depicted in Fig. [T]), L(a\M{) is the so 
called log-likelihood functional 

n 

L(a\Mi)=J2^og Pj (Mi\a), (12) 

i=i 

and Pj(Mi\a) denotes the conditional probability for the j-th run. We observe that Eq. (fTTj) is concave since the 
logarithm of a linear function is a concave function and the summation of concave functions is a concave function. 

To solve the ML problem, we employ the iterative numerical method described in First, we generate a large 
amount of data distributed according to Eq. ([8]), for some fixed value of a and j3. Then, we fix some order zero 
approximation ao for the estimator qml- Then, the first order correction is given by 



dL(a|Mi) 




da 






a— ao 



- - FM <13) 
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where F(a) is the Fisher information 



which in our case is given by 



F[a) = S ( mJ RmTh ' (15) 



Mi=0 



The Fisher information measures the amount of information that the random variable Mi carries about the unknown 
parameter a on which the likelihood function depends. So, the estimator to first order is ao + a\, and the procedure 
can be iterated with this value as order zero approximation to obtain higher order corrections. Obviously, the result 
is independent of the initial value ao- For f3 not too big (say < j3 < j), the algorithm converges in a few steps (say, 
less than 10). 

The variance on the ML estimator of a parameter satisfies the Cramer-Rao bound [HI 

o 2 {a Mh ) > — \— (16) 
nt (a) 

The bound in Eq. (|16p is saturated if the number of data is large enough and the parameter is mono-dimensional (as 
in the present case). We numerically estimated the variance of the estimator Qml in Eq. (jlip by dividing the data 
into blocks, finding the estimator aj for each block, and then calculating the variance of such estimators, namely 

v 2 (aML) = ^(ai - a ml) 2 - (17) 

i 

In Fig. [4] we verified that the variance numerically saturates the bound in Eq. (|16[) . 




FIG. 4: (Color online) Variance a 2 (aML) versus the number n of outcomes (runs), for parameter a — 0.75, measurement noise 
j3 = 0.25, and number of POVMs M = 10. The solid line represents the Cramer-Rao bound in Eq. (I16|l . 

Fig. [5] shows that the variance decreases as the number of POVMs used in parallel increases, and upper and lower 
bounds for variance. To find the upper bound consider the function 

/(Mi) - , (18) 

Notice that /(Mi) is an unbiased estimator for the parameter a, since one has 



M 

(f(Mi))= f(M 1 )p{M 1 \a) = a. 

Mi=0 



(19) 
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FIG. 5: (Color online) Variance a 2 {aML) versus the number M of POVMs, for parameter a — 0.75, measurement noise f3 = 0.25, 
and number of runs n = 2000. The upper and lower bound correspond to Eqs. (|21[) and (|22[) . respectively. 



The second moment is given by 



A/ 



(/(Afx) 2 ) = f{M 1 fp{M 1 \a) 

Mi=0 

(l-2/3) 2 Af 

Thus, an upper bound for the variance on the parameter a is 

a{a ML )<[a-a + (1 _ ^ ) ~ 



(20) 



The lower bound for the variance is 



<y 2 {a M L) > 



1 



nF(a)\f3 =0 ,M= 



(21) 



(22) 



where the right-hand side of Eq. (|2"2"|) corresponds to the use of the ideal detector on the original state. 

The computed variance shows a dependence on the parameter a similar that in Eq. (|2ip . decreasing as 1/Af for 
a = or a = 1. The variance saturates the lower bound in Eq. (1221) . so the the estimator of the parameter a is 
optimal. 



IV. MAXIMIZATION OF MUTUAL INFORMATION 

We consider now the mutual information as the figure of merit in the measurement purification scheme. We consider 
a qubit parametrized as 



cos-|0) +e # sm-|l) 



The probability in Eq. can be written as 

( M 



p{M x \9) = 



[(1 - (3) Mo (3 Ah - (1 - /3) Ah f3 Mo ] cos 2 ~ + (1 - /3) Ml /3 Mo 



(23) 



(24) 



independent of • 
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In the following we suppose that the prior probability p(9, <fi) of having the input state in Eq. (|23|) is uniform, so 
that the mutual information I{Mi : 9, (f>) between random variables 9 and </> and random variable M\ is given by 

/(Mr : 6) := ~ C dBsmO £ pfM^) log ( i rJ^fLim ) ' ^ 
The integral in the denominator gives 

£ de S mBp(Mi\B) = (^)((l-/3) Mo /5 Ml + /3 Mo ), (26) 

and a lengthy analytical form for Eq. ()25[) is provided in the Appendix of the paper. 

The mutual information I{M\ : 9) saturates the bound I{M\ : 9) < I{M\ : 9)\m=i,/3=o — 0.279 bit for increasing 
M. Notice that the mutual information does not converge to 1 bit, since a continuous "alphabet" of states is allowed. 
The mutual information /(Mi : 9) saturates almost exponentially versus M , as shown by Fig. [6l This means that we 




FIG. 6: (Color online) Function 
The mutual information 7(Mi : 9 



-log, 1 HM 11 0)_ 

is given by Eq. (|A1|I 



j versus the number M of POVMs, for measurement noise /? = 0.25. 



are recasting the use of many noisy detectors to an effective use of a single ideal detector. 

Let us consider in more detail the simplified case in which the only allowed input state are the up (9 = tt) and 
down (9 = 0) eigenstates of a z . This simplification leads to two advantages: a much more tractable analytical form 
for the mutual information l2(Mi : 9), and the possibility to make a comparison with classical post-processing based 
on majority voting. The mutual information is given by 



M 



M 
Mi=0 v 



(1 _ p)MipMo bg 



2(1 -p) M ^ 



Mi oM Q 



(1 - /3) M °/3 M i + (1 - /3) M i/3 M ° 



(27) 



Eq. ([2~7T) behaves as expected for the ideal POVM case (i.e. (3 — 0), where ^(Mx : 6) = 1, and for completely 



isotropic POVM case (i.e. j3 = |) where 



0. Finally, we investigate the optimal classical post-processing 



to be applied on the M outcomes of the parallel noisy POVMs to maximize the mutual information. We simply argue 
that majority voting is close to the optimal post-processing, as is shown in plot Fig. [71 The gap between the binary 
mutual information li{M\ : 9) and that obtained with majority- voting strategy could be explained by the fact that 
in general a number of POVM elemen ts g reater than the cardinality of the input alphabet can optimize the mutual 
information. In fact, Davies' theorem [16| puts an upper bound of d 2 on the number of POVM elements to optimize 
the mutual information for an alphabet of d linear independent pure states (see also [ll|h 

The case of two-states alphabet can be easily generalized to an alphabet of d orthogonal state {|j}}, with j — 
1,2, .. .d, and noisy POVM elements P[ = + =-jjrI- The conditional probability of the outcomes of M noisy 

measurements on M copies of \j) is simply the multinomial 



p(M 1 ,M 2 , ■ ■ ■ , M d -i\j) 



Ml 



[(d-l)a + l] M i{l-a) 



\M-Mi 



M 1 \M 2 \---M d 



d M 



(28) 
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FIG. 7: (Color online) Function — log 2 [l — Ii{M\ : 8)] versus the number M of POVMs, for measurement noise /3 = 0.25. 



where Mi is the number of outcomes I in the string of M outcomes, and Ma = M — Ylj=i Mj. The conditional proba- 
bility allows one to evaluate the mutual information, and for increasing number of clones M, the noisy measurements 
are purified. In Fig. [8] we show the purification effect for an alphabet of four orthogonal equiprobable states, and two 
different values of noise. As expected, for increasing value of M the mutual information approaches two bits. 




FIG. 8: (Color online) Mutual information for an alphabet of four equiprobable orthogonal states versus number of purifying 
copies, for measurement noise a = 0.8 (upper), and a — 0.4 (lower). 



V. INEFFICIENT PHOTODETECTION 



In the rest of the paper we consider observables with unbounded spectrum, for which the concept of amplification 
applies, and we review the scheme of H. P. Yuen [l[ for improving noisy photodetectors, homodyne and heterodyne 
detectors. In the original proposal of Ref. [l[ the signal-to-noise ratio improvement was studied for noisy measurements 
with preamplification assistance. By reviewing the results here, we explicitly consider the effect of amplification as a 
purification of the noisy POVMs, and hence of the outcome probability distributions. 

Light is revealed by exploiting its interaction with atoms/molecules or electrons in a solid, and, essentially, each 
photon ionizes a single atom or promotes an electron to a conduction band, and the resulting charge is then amplified 
to produce a measurable pulse. In practice, however, available photodetectors are not ideally counting all photons, 
and their performances is limited by a non-unit quantum efficiency rj, namely only a fraction r\ of the incoming photons 
lead to an electric signal, and ultimately to a count: some photons are either reflected from the surface of the detector, 
or are absorbed without being transformed into electric pulses. Let us consider a light beam entering a photodetector 
of quantum efficiency rj, i.e. a detector that transforms just a fraction 77 of the incoming light pulse into electric 
signal. We will focus our attention to the case of the radiation field excited in a stationary state of a single mode at 
frequency w. Then, the Poissonian process of counting gives the following probability Prj(m) of revealing m photons 



Pr ,(m) = Tr L (r?a f m expC-^a):! , (29) 
ml 
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where p represents the quantum state of light, and : : denotes the normal ordering of field operators. 
Using the identities 



:{a}a) n : = (a^) n a n = a)a{a)a - 1) . . . {a* a - n + 1) , (30) 



. e -xa"a : = ^ _Z£L( a iy a l = (I _ a .)a t a j ( 31 ) 
Z=0 



one obtains 



where 



n— ?n 



p,(m) = £ p„„ ( " ) v m {l ~ vT~ m , (32) 



Pnn = (n\p\n) = p n= i(n) . (33) 

Hence, for unit quantum efficiency a photodetector measures the photon number distribution of the state, whereas 
for non unit quantum efficiency the output distribution of counts is given by a Bernoulli convolution of the ideal 
distribution. 

The outcome distribution in Eq. (|32[) can be equivalently described by means of a simple model in which the 
realistic photodetector is replaced with an ideal photodetector preceded by a beam splitter of transmissivity r = rj. 
The reflected mode is absorbed, whereas the transmitted mode is photodetected with unit quantum efficiency. In 
order to obtain the probability of measuring m clicks, notice that, apart from trivial phase changes, a beam splitter 
of transmissivity r affects the unitary transformation of fields 



where all field modes are considered at the same frequency. Hence, the output mode c hitting the detector is given 
by the linear combination 



c = \fra — \f\ — rb , (35) 

and the probability of counts reads 

p T (m) = Tt[U T (p®\0)(0\)l4\m)(m\®l] 

(l-T) n - rn T m . (36) 



^ p ™ ( m ) 



Equation (|32p is then reproduced for r = 77. Wc conclude that a photodetector of quantum efficiency rj is equivalent 
to a perfect photodetector preceded by a beam splitter of transmissivity 77 which accounts for the overall losses of the 
detection process. According to Eq. (|3"2"1) . the POVM describing the inefficient photodetector can be written as 



n„(m)= ^ \ v m (l-r,y ta - m , (37) 

such that Pri(m) = Tr[pIL n (m)]. The random variable m, suitably rescaled by 77, provides an estimator of the average 
photon number (a^a) = Tr[pa'a], since one has 



m 

m— 

In order to evaluate the second moment of the probability, one uses the identity 



V — n„(m) =a f a . (38) 



00 / m \ 2 1 _ 

V( — I UJm) = (a f a) 2 + -a) a , (39) 

V * J 7 ? 

and hence the inefficient measurement is affected by the added noise — -{a? a), with respect to the ideal intrinsic noise 
A(at a ) 2 = ((aM 2 ) - {aU) 2 . 
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In the following we show that and ideal photon-number amplifier can arbitrarily reduce the added noise of the 
inefficient measurement for increasing gain. The ideal photon-number amplification map is given by [l8l - |20j 

a f a — >V f a)aV = go) a , (40) 

where g is an integer, and V is the isometry 



V = J2\gn){n\. (41) 



n=0 

The preamplihed POVM is simply given by 



II^(m) = ( 9a ^ a ) if\l - v r' a ~ m . (42) 
The estimator of the average photon number is now m/(gr]), and the second moment is given by 



E 

m=0 



—) 2 I#) ( m ) = ( a t a) 2 + ^—^a)a . (43) 
9Vj 9V 



Clearly, for g — > oo, the added noise is completely removed for any value of the quantum efficiency rj. 

We notice that the ideal photon-number amplifier is so effective that indeed even a preamplified heterodyne detection 



provides the ideal photon number distribution for increasing gain, as shown in Ref. [21 



VI. INEFFICIENT CONTINUOUS VARIABLE MEASUREMENTS 
A. Homodyne detection 

The balanced homodyne detector provides the measurement of the quadrature of the field 

X, = + ae (44) 

It was proposed by Yuen and Chan :22j], and subsequently experimentally demonstrated by Abbas, Chan and Yee 
[23j |. The signal mode a interferes with a strong laser beam mode 6 in a balanced 50/50 beam splitter. The mode b 
is the so-called the local oscillator (LO) mode of the detector. It operates at the same frequency of a, and is excited 
by the laser in a strong coherent state \z). Since in all experiments that use homodyne detectors the signal and the 
LO beams are generated by a common source, we assume that they have a fixed phase relation. In this case the LO 
phase provides a reference for the quadrature measurement, namely we identify the phase of the LO with the phase 
difference between the two modes. By tuning (p — argz we can measure the quadrature X v at arbitrary phase. 

Behind the beam splitter, the two modes are detected by two identical photodetectors (usually linear avalanche 
photodiodes), and finally the difference of photocurrents at zero frequency is electronically processed. In the strong-LO 
limit \z\ — > 00, the homodyne detector is described by the POVM 



dX 
2tt 



U ( x ) = I exp[U(X v - *)] = \x) vtp {x\ , (45) 



oc 



namely the projector on the eigenstate of the quadrature X v with eigenvalue x. In conclusion, the balanced homodyne 
detector achieves the ideal measurement of the quadrature X v in the strong LO limit. In this limit, the probability 
distribution of the output photocurrent approaches exactly the probability distribution p(x,<p) = v (x\p\x} v of the 
quadrature X v , and this for any state p of the signal mode a. 

It is easy to take into account non-unit quantum efficiency at detectors. The POVM is obtained by replacing 



X - r X. - \j-^-( u v + v v) ( 46 ) 
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in Eq. (|45|) . with w v — (w'e l<p + we lip )/2, where w — u,v denotes the vacuum modes of the two inefficient 
photodetectors. By tracing the vacuum modes u and v, one obtains 



+ 00 



— OO 

1 



dX 
2^ 

dX 
2^ 



(47) 



: exp 



(x - X v ) 
2A2 



2 1 



+ OO 1 / _ f\2 

dx'e X> \x')^{x'\ 



where 



A,, 



4?7 



(48) 



Thus the noisy POVM, and in turn the probabi lity distribut ion of the output photocurrent, are just the Gaussian 
convolution of the ideal ones with r.m.s. A r; = \/(l — rj)/(4r]). 

In the following we show that the added noise of the inefficient homodyne detector can be removed by amplifying 
the signal by means of a phase-sensitive amplifier. This amplifier is described by the squeezing operator 



S(H) = cxp 



(49) 



and performs the mode transformation 



&(£)aS(£) = (cosh|^|)a+ ||(sinh|C|)a 



For ^ = re 2ltp , with r > 0, one has 



(50) 



(51) 



Hence, the effective POVM obtained by preprocessing H v (x) in Eq. (|4"5|) with the phase-sensitive amplification of X v 
is given by 



n< r »(i) = S\re 2 ^)Ii v {x)S{re 2 ^) 



27TA2 



: exp 



(x - e r X v f 

2A 2 

v 



(52) 



Now, in order to obtain an unbiased measurement of X v , it is enough to rescale the outcome by e r . On the other hand, 
the added noise with respect to the ideal measurement X v becomes equal to e~ 2r A 2 , which can be made arbitrary 
small for increasing value of the squeezing parameter r. 



B. Heterodyne detection 

Heterodyne detection allows to perform the joint measurement of two conjugated quadratures of the field [24l [25|. 

A strong local oscillator at frequency w in a coherent state \a) hits a beam splitter with transmissivity r —> 1, and 
with the coherent amplitude a such that 7 = \a\ \J r(l — t) is kept constant. If the output photocurrent is sampled 
at the intermediate frequency ujif, just the field modes a and b at frequency w ± u)if are selected by the detector. 
Modes a and b are usually referred to as signal band and image band modes, respectively. In the strong LO limit, 
upon tracing the LO mode, the output photocurrent I(ujif) rescaled by 7 is equivalent to the complex operator 

Z=!^lA = a-b\ (53) 
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where the arbitrary phases of modes have been suitably chosen. The heterodyne photocurrent Z is a normal operator, 
equivalent to a couple of commuting self-adjoint operators 

Z = RcZ + i\mZ , [Z, Z 1 ] = [RcZ, ImZ] = . (54) 

The POVM of the detector is then given by the orthogonal (in Dirac sense) eigenvectors of Z. 

In conventional heterodyne detection the image band mode is in the vacuum state, and one is just interested in 
measuring the field mode a. In this case the POVM H(z) is obtained upon tracing on mode b, and one has the 
customary projectors on coherent states 

U(z) = -\z){z\, (55) 

with z E C. The coherent-state POVM provides the optimal joint measurement of conjugated quadratures of the field 
5]. For a state p, the expectation value of any quadrature X v is obtained as 

(XJ = Tr[pX v } = / —Re(ae-^)(a\p\a) . (56) 
Jc n 

The price to pay for jointly measuring non-commuting observables is an additional noise. The r.m.s. fluctuation is 
evaluated as follows 

f ^[Re(a e -^)] 2 (a|p|a) - (V^) 2 = (AV 2 ) + - , (57) 

where (AX 2 ) is the intrinsic noise, and the additional term is usually referred to as "the additional 3dB noise due to 

the joint measure" [HHH- 

The effect of non-unit quantum efficiency can be taken into account in analogous way as in Sec. I VI Al for homodyne 
detection. The coherent-state POVM is replaced with the convolution 



d 2 z' -^^\z'){z'\ 

7TA2 



/ ZXI e ^T^^, (58) 



v 



where A 2 = (1 — ?/)/??. 

In the following we show that inefficient heterodyne detection can be purified by phase-insensitive amplification. 
Phase-insensitive amplification with (power) gain G amplifies the coherent amplitude of coherent states by vG, at 
the expense of addition thermal photons n = G — 1. Differently from phase-insensitive amplification, the physical 
process is not unitary, but described by a completely positive map £q. Here, we just need the Heisenberg evolution 
of the projector on coherent states, which is simply given by the rescaling [29] 

^(\a)(a\) = ^\G- 1 / 2 a)(G- 1 / 2 a\. (59) 



It follows that under phase-insensitive preamplification the noisy heterodyne POVM (|58|) is replaced with 

(60) 

Upon rescaling z — > G~ 1//2 z, one obtains 

Gd 2 z' -°^^*\ z '){z'\ 



7TA2 



■e 



\z){z\ 

~G— >CXD 



(61) 



namely the noise due to quantum efficiency can be arbitrarily reduced for increasing value of the gain G. The 
effectiveness of preamplification in purifying heterodyne detection is more unexpected than the case of homodyne 
detection, since phase-insensitive amplification is not a unitary process. 
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VII. CONCLUSION 

In this paper, we addressed the problem of optimizing a quantum channel acting before many parallel uses of a 
noisy POVM in order to purify the measurements, namely to achieve an effective measurement that is less noisy than 
the original ones. We first considered the purification of a z noisy measurements on qubits, by choosing the orthogonal 
cloning channel as a purification map. We found the maximum-likelihood estimator for o~ z , whose variance shows 
an almost-exponential decay of versus the number of POVMs. We also worked out an analytic form for the mutual 
information between the state parameter and the outcomes of the POVMs, and here also an almost-exponential 
improvement versus the number of POVMs has been found. We proved that naive majority voting is not the optimal 
classical post-processing, since the maximum-likelihood approach gives a better estimator. For photodetection and 
continuous variable measurements as homodyne and heterodyne detection, the measurement purification can be 
achieved by preamplification, as early pointed out by H. P. Yuen [l|. 

We think that the relevant problem of purifying noisy quantum measurements will have a significant impact on the 
quantum information technology, in this same way as the decoherence problem. 
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Appendix A: Derivation of the mutual information in Eq. (|25[) 

We provide here an analytic form for Eq. (|25|) . obtained substituting Eq. (jHJ for p{M\\6). One obtains 

Mi=0 v ' 

where 

ci = (1 - /3) 2M /? 4Ml , c 2 = (1 - /3) 4Ml /3 2M , (A2) 

and 

A = - 161og((l - pf M ^ M + (1 - p) M 0* M *) 



+ 81og((l - /3) 2M /3 4Ml ) + 16 - 8-i-. 

ln(2) 



(A3) 



B = + 161og((l - /3) 2Ml /3 M + (1 - f3) M f3 2Ml ) 

- 8 log((l - /3rV M ) - 16 + 8-i-, (A4) 

ln(2) 

C = (1 - p^h _ (1 _ ^3Mi pM+M, (A5) 
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